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Abstract. Within the variational approach of Bose liquids we analyze the g.s. energy of charge 
neutral alpha matter at T=0. As a prerequisite for such calculation we take from the literature or 
propose new a — a potentials that are particularly suitable for this task, i.e. posses a repulsive 
core and/or reproduce the low energy scattering data and the resonance properties of the a — a 
system. The alpha matter EOS is then obtained with the HNC method using Pandharipande-Bethe 
correlation derived variationally in the lowest order expansion of the energy functional or a simple 
gaussian function with a healing range determined by the normalization of the radial distribution 
function in the lowest order. We show that saturation is achieved only via repulsive and shallow 
potentials that are not consistent with the scattering and resonance constraints. 
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INTRODUCTION 

The pioneering alpha-matter calculations reported in [1] made use of a — a potentials 
characterized by a strong or even infinite repulsive component. These interactions were 
constructed to fit the elastic- scattering phase shifts deduced from experiment. The Equa- 
tion of State (EOS) depends strongly on the shape and strength of these potentials and 
it is usually derived within the frame of the paired-phonon analysis (PPA) or the HNC/n 
method. The EOS calculated with hard-core potentials, trivially saturates at densities 
and energies close to the nuclear matter saturation point (p a w 0.04 fm~ 3 , E/N a ~ -11- 
16 MeV). On the other hand, for a typical soft-core potential such as the one proposed 
by Ali and Bodmer (AB) [2] the alpha matter almost fails to saturate [1]. In fact a deep 
minimum in a very soft EOS at a high density is predicted with AB potential. However 
at such high densities the alpha-condensate is almost completely depleted due to Pauli 
blocking. Somehow this disappointing result is conflicting with what one would expect 
based on the manifestation of alpha clustering in real nuclei. The clusterization of alpha 
particles on the surface of nuclei at densities around half the central nuclear density, 
as revealed by a-decay, a-transfer reactions or the putative dilute three-alphas conden- 
sate in the Hoyle state of 12 C, are pointing to a high stability of alpha matter at lower 
densities. 

We recently addressed the problem of alpha matter saturation and concluded that in 
order to avoid the collapse of the EOS very shallow potentials with a strong repulsive 
core are necessary [3]. It would be then of interest to extent this study to other potentials, 
especially those in agreement to elastic phase-shift data or the first resonant states in the 



short lived Be. It is therefore timely to gain a better understanding of the role played by 
the a — a potentials in a many-body approach to a-matter. 

ALPHA-ALPHA POTENTIALS 

The popular 5-state Ali-Bodmer a — a potential consists of a short-ranged (1.43 fm) 
repulsive part and a long-ranged attractive part (2.50 fm) [2], 

Vaa(r) = 475exp [-(0.7r) 2 ] - 130exp [-(0.475r) 2 ] . (1) 

This potential reproduce the a — a elastic scattering phase-shifts at low energies. The 
same phase-shift properties are shared by the deep Buck et al. (BFW) potential [4] , 

Vaa(r) = Vbexp[-Gur) 2 ], (2) 

where V =-122.6225 MeV, ji = 0.469 fm" 1 . In addition the BFW potential reproduces 
the energy and the width of the first + resonant state in 8 Be. 

A microscopic approach to the a — a potential is provided by the double-folding 
method [5] where the interaction between two alpha ions is calculated as an overlap 
of the single particle densities of the interacting objects smeared by the local two-body 
potential v nn , 

V aa (r) = J dr x J dr 2 pa(ri)pa(r 2 )v nn (p ,r - r\ + r 2 ) (3) 

The effective n — n interaction v nn is taken to be dependent on the density p of the 
nuclear matter where the two nucleons are embedded. It should also consist of a density 
independent finite-range part with preferably two ranges such that a potential similar to 
the Ali-Bodmer is obtained at least for the direct component. A choice satisfying these 
requirement is provided by the Gogny [6] interaction. Very recently we proposed a bare 
a — a interaction based on the double-folding method at energies around the barrier, 
using realistic densities of the a -particle and the density dependent Gogny nucleon- 
nucleon effective interaction [3]. The bare a — a potential consisted only in the direct 
term. In the present work we add the knock-on nucleon exchange term including recoil 
corrections and the non-local kernels are localized in the lowest order of the Perey-Saxon 
approximation at energies around the barrier [7]. Since all three parametrizations of the 
Gogny force (Dl [8], D1S [9], DIN [10]) are dominated by exchange the resulted folded 
potentials are very deep. It is a matter of evidence that deep potentials such as BFW and 
Gogny are never saturating alpha matter because of the large amount of attraction that 
will contribute to the total energy. 

For our alpha matter investigations let us first consider the direct part in the double- 
folding potential (3). The motivation behind such a simplification was recently discussed 
in connection with the necessity of accounting for the incompressibility of nuclear matter 
in cold clustering processes [11] and extreme sub-barrier fusion [12]. In this framework 
a double-folding repulsive potential with a zero-range interaction is added to the direct 
and exchange potential such that the energy cost for overlapping two pieces of nuclear 



matter are payed-off. The strength of this repulsive 8 -like potential is in a simplified 
picture proportional to the nuclear incompressibility at the corresponding density of total 
overlap. 

We assume a Gaussian nuclear matter distribution of the a -particle, 
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where the length parameter b is determined from the root mean square radius (rms) 
1.58±0.002 fm extracted via a Glauber analysis of experimental interaction cross sec- 
tions [13]. 

Using the analytical expression of the direct component n — n force in the Gogny 
parametrization, as listed in [3], and inserting the Gaussian density distribution (4) in 
the double folding integral (3) the direct part of the a — a interactions reads, 
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Above, Wi,Bj,Hi,Mi are strength of the i = 1, 2-th finite-range term of the Gogny 
parametrization whereas jU ; the corresponding ranges, while the parameter y charac- 
terizes the density dependence of the force. In the left panel of Fig.l we represent the 
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FIGURE 1. Comparison of the Gogny Dl and DIN direct a — a potentials with the AB potential. 

AB, Dl and DIN potentials on a magnified scale around the minimum. The two Gogny 
potentials display pockets that are shallower and shifted to larger radii compared to the 
Ali-Bodmer potential. The minima of these potentials are close to the classically touch- 
ing configuration of the a — a system. 

The derivation of the exchange part of the potential will be given elsewhere. Here we 
only reproduce the closed expression obtained after the localization via the Perey-Saxon 



method of the non-local kernel, 
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Above we distinguish between a sub-barrier branch of the potential (K 2 <0) and an 
over-barrier one (K 2 >0). Note that in the exchange kernel we used an alpha single- 
particle density matrix constructed from the Os harmonic-oscillator orbitals. Only a few 
iterations are needed to obtain convergence in the localization procedure. In the left 
panel of Fig. 2 the full (direct+exchange) a — a potential for three parametrizations of 
the Gogny interaction is displayed at zero relative momentum (X rel = 0) and compared 
to the BFW potential. 



Supersymmetric Potentials 

In a previous work we concluded that a shallow potential (in any case shallower than 
the traditional AB) is needed in order to obtain the saturation alpha matter at physically 
reasonable densities [3] . On the other hand, a shallow potential results naturally from the 
requirement that 8 Be has no bound states. An alternative way to the Ali-Bodmer phe- 
nomenological ansatz to construct a shallow potential was proposed by Baye [14]. A se- 
quence of supersymmetric transformations, which removes the Pauli-forbidden states in 
the sense of orthogonality condition model (OCM), is applied to the original Schrodinger 
equation in such a way that the phase-shifts corresponding to the transformed Hamilto- 
nian coincide to the phase-shifts of the original Hamiltonian. In this manner one obtains 
a new shallow potential that is exactly phase-equivalent (modulo 2k) to the original 
deep potential. In the particular case of a — a system the phase-equivalent potential is 
obtained from the suppression of the two lowest (/=0) bound forbidden states. A com- 
pact formula provides the needed shallow potential 

2 2 

V*2(r) = V£2(r) " -^logdct¥ 2 \r) (9) 



where the elements of the 2x2 matrix read 

(^ 2) )ij(r) = f r u i (r')u J (r')dr f (ij = 0s,ls) (10) 

and Vaa is any one of the deep potentials discussed earlier. Above we denote by uq s xs the 
lowest two bound state wave functions of the deep nuclear+Coulomb (/=0) potential. We 
obtain these states using the Runge-Kutta-Nystrom (RKN) numerical integration [15]. 

The suppersymmetric shallow partners of the BFW potential were determined already 
by Baye [14] and more recently in a semi-classical approach in [16] for 1=0, 2 and 4. 
The first two bound states in this potential are located according to this reference at - 
72.8 MeV and -25.9 MeV. With the RKN method we obtain the values -72.62 MeV and 
-25.61 MeV. The minimum of the shallow BFW-SUSY potential is located at 2.85 fm 
and has a depth of -7.56 MeV and is similar to the Baye estimation and about 1 MeV 
shallower than the semiclassical estimation of Horiuchi et al. [16]. 

The transformed potential displays a 20/r 2 singularity at the origin. In practical 
calculations we use an interpolation with gaussians that is non-singular at the origin. 

To obtain the SUSY partner of the Gogny potentials we constrain first the orig- 
inal deep potential to reproduce the properties of the first resonant state in 8 Be, i.e. 
E KS =92.12 ± 0.05 keV, ro=6.8± 1.7 eV [17]. The renormalized Gogny potentials are 
next subjected to a sequence of two supersymmetric transformations in order to elimi- 
nate the Pauli-forbidden states. In the right panel of Fig. 2 we display the supersymmetric 
partners of the BFW and the three Gogny interactions (including the Coulomb compo- 
nent). We make the remark that the deep potentials BFW, Dl and DIS(DIN) have very 
different depths but upon renormalization and supersymmetric transformations all four 
potentials have similar depths. Thus, the physical contraint imposed by the + resonance 
in 8 Be and the removal of the forbidden states, lead to an almost unique potential for the 
a — a system. 

EOS OF ALPHA MATTER 

We derive the EOS of alpha matter using the variational approach, described in [3]. In 
Fig.3 we compare the EOS for AB potential with BFW-SUSY and Dl-SUSY and we 
conclude that for this last two cases the EOS curves are collapsing, while the EOS for 
the AB interaction saturates at too large densities. Thus shallow potentials constrained to 
reproduce the phase-shift elastic data and/or the first resonant state in 8 Be do not saturate 
alpha matter in the HNC/0 approximation. The reason for this disappointing result is that 
the two-body correlation functions obtained with the Pandharipande-Bethe prescription 
display large overshootings near the healing distance with increasing alpha matter den- 
sity. The density dependence of the healing distance d does not bring sufficient density 
dependence in order to ensure saturation. 

In order to demonstrate this effect we consider a simple phenomenological two-body 
Gaussian correlation function [18], 
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FIGURE 2. Bare direct+exchange a — a potentials for the BFW and the three n — n Gogny interactions 
used in the text (left panel) and their supersymmetric partners (right panel). 
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FIGURE 3. Comparison between the EOS predicted with HNC/0 for AB, BFW-SUSY and Dl-SUSY. 



Evidently this function has no overshooting. The parameter /3 is determined from the 
normalization condition [19], 

POO 

Anp / drr 2 (f 2 (r) -1) = -1 (12) 
Jo 

where p is the density of alpha matter. We perform a selfconsistent HNC/4 calculation 
using these correlation functions and the bare AB potential. Both Jackson- Feenberg 
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FIGURE 4. Equation of state for alpha matter in the HNC/4 approximation. The bare AB potential is 
scaled with a factor N in order to see the variation of the saturation point as a function of the potential 
strength. 



(Ejf) and Pandharipande-Bethe (Eps) are displayed in Fig. 4 for different renormaliza- 
tions N of the AB potential. The intrinsic binding energy of the a particle (B a ~ 7 MeV) 
is not included. The PB energy is calculated using the Kirkwood approximation for the 
three-body radial distribution. The values obtained in this approximation are quite simi- 
lar with those obtained with JF approximation since the contribution from the three-body 
term is negligibly small. The saturation density does not change with the renormaliza- 
tion of the potential since this is determined exclusively by the density dependence of 
the range parameter /3 ~ p 1 / 3 (Eq.12). Similar results have been obtained for the equa- 



tion of state using a 4 order cluster expansion and the HNC/0 method. Regardless of 
the normalization factor the EOS saturation point is very close to the saturation point of 
normal nuclear matter, i.e. Pnm = 4p« « 0.16 fm~ 3 . 

In conclusion, shallow phenomenological and semimicroscopic a — a potentials 
which reproduce low energy a — a scattering phase shifts and/or the resonance in 8 Be 
lead to a collapsing alpha matter equation of state if two-body correlation functions are 
obtained with the Pandharipande-Bethe prescription. In this prescription, the assumption 
of a finite weakly density dependent healing distance leads to overshooted two-body cor- 
relation functions and enhances the importance of the long attractive tail of the potential. 
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